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Abstract 

On the basis of Bethe ansatz solution of one dimensional SU(4) Hubbard- 
like model, we study its thermodynamics properties by means of Yang- Yang 
thermodynamics Bethe ansatz. The Lande g factor is taken into account 
so as to describe electrons with orbital degeneracy. The free energy at low 
temperature is given and the specific heat both in strong coupling and weak 
coupling limits are obtained. 

I. INTRODUCTION 

The study of integrable model has a long and rich history in condensed matter physics 
beginning with Bethe's solution of the one-dimensional Heisenberg Model and extending to 
nowadays when a variety of soluble models provide the paradigms that enrich much of our 
physical intuition. Recently there has been much interests in the studies on the Sci electrons 



in transitional metal oxideS because of the existence of orbital degree of freedonJli in addition 
to the spin one. Since the SU(4) symmetry was pointed oulJl, there has been various studies, 
such as critical properties in photoemission spectra for the one-dimensional Mott insulator 
with orbital degeneracyi. One dimensional solutions of the model in insulating limiti as well 
as the Hubbard-like mo deli were presented. There were discussions of the thermodynamics 
based on the SU(N) generalized Bethe ansatz equation^, but the physical meaning of the 
external field is ambiguous there because it is not associated with spin and orbital directly. 

In present paper, we study the low-temperature thermodynamics of SU(4) Hubbard- 
like model for electrons with 2-fold orbital degeneracy at low concentration. The thermal 
equilibrium is discussed exactly on the basis of the known Bethe-ansatz equations by tak- 
ing account of the Lande g factor. Next section we first briefly exhibit the Bethe-ansatz 
equation under consideration, then discuss the thermodynamics of the model by considering 



string hypothesis. In Sec. |T|, we calculate the formal expressions of the thermodynamics 
quantities, such as free energy etc.. In Sec. 0, the case of low temperature is discussed 
extensively. 



II. THE MODEL AND ITS SPECTRUM IN THERMODYNAMICS LIMIT 

The model Hamiltonian reads 

^ = E nCtaC^+l,a + Ct+,^aCi,a)V + f/ ^ ^.a^i,^' (l) 
i,a i,a<a' 

where i = 1,2, identify the lattice site, and a = 1,2,. ..,4 labels the four states of 

spin and orbitali. The V projects the Hilbert space onto the sector that the sites are only 
occupied by at most two electrons. The internal degree of freedom (|I]) is specified to spin 
and orbital in present model. The Bethe-ansatz equation for the spectrum was suggestedi'0 



as 



,^k,L ^ -Q sin kj -Xb + iv 



^ sin kj — Xb — ir]^ 
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SL Xg - sin ki + if] _ _ Xg - Xb + i2r] ^ /ic - Ag + ii] 
1=1 ^'^ ~ ~ 6=1 Aa - Afe - i2ri f J fi^ - Xg - ir]' 

-pi /ig - Ab + i?7 _ _ yr fJ'a- fJ'c + i2r] |C' - /ig + ir] 

6=1 f^a- Xb- it] fJ'a- f^c- i^T] Vd ^ - 

bJ{ Va- fJ-b- ir) ]^c- i^T] ' 

where t] = U/At. Unlike the conventional SU(2) Hubbard model the SU(4) generalization 
(@) of Lieb-Wu solution! is valid at low temperature and low concentration00. 

For ground state (i.e., at zero temperature), the fc, A,/i, are real roots of the Bethe 
ansatz equation (0). However, for the excited state, they can be complex roots. We will 
not take account of the complex roots in the charge space k for repulsive interactionEl. The 
complex root of A, fi, v are always form a "bound state" with several same kind of rapidity 
which arise from the consistency of both hand side of the Bethe ansatz equations. The 
complex roots are defined as 

A:^' = a: + + 1 - 2j>r7, j = l,2...,n, 
t/,'"^- = + (m + 1 - 2j)^r7, j = l,2...,m, 

Vi' = v\^{l^\-m"n. J = l,2...,/, (3) 

up to the order 0(e~^) which vanishes in the thermodynamic limit. Substituting those 
strings into Eq. (H) and taking the logarithm of it, we get 



27rL = k,L + 2 V tan-i f ^" 



nr] 

I \n 



2nJ: = 2 5: tan-i _ 2 J] tan'^ (^M^) " ^ E BnU tan"^ (^^^) , 



n \ I Ti I in 

2nK: = 2 E B^u tan-^ f^^^) - 2 E ^mt tan'^ fl^lL^\ _ 2 ^ B^u tan"^ (^^ ~ 

tT] 17] 
n 1,1 i,n ,,l 



clt 



27rg: = 2 E 5n^* tan-i (^^^^) - 2 E ^n^* tan'^ • (4) 

where 
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A 



nmp 



Bnlt — ' 



1, for p = m + n, |m — n| (7^ 0), 

2, for p = n + m — 2, n + m — 4, • • • , |n — m| + 2, 

0, otherwise, 

1, for t = n + / — 1, n + / — 3, |n — /| + 1, 
0, otherwise. 

Noe we consider the question in the thermodynamic hmit N, L, M, M', M" — > 00 with a 
fixed concentration D — N/ L by introducing the distribution of roots and holes for /c, A, v 
respectively: 



1 dKn{tj) 
L dfi 

L dv 



(5) 



We obtain the following coupled integral equations. 

p + p'* = 7^ + Vcos/c / Kn{sva.k - X)an{X)d\, 
a„ + = / Kn{\ - smk)p{k)dk - ^A^mi f Ki{X - X')am{X')dX' + Y.Bnit f Kt{X - ii)uJi{ii)dii, 

■' ml ■' It ■' 

cUn + uJn = ^Bnit I Kt{fi - X)ai{X)dX -J2^nH Kt{n - n')uJi{n')dij! + Y,Bnit I Ktifi - u)n{u)du, 
It •' It It 

+ r!:^Yl Bnlt I Kt{v - iJi)ui{ix)dix - ^ A,,n [ - u')Ti(u')du' . (6) 
It ^ It ^ 

where Kn{x) — T:~^nr]/ {ri^rf + x^), and the integration limits are defined by 

— = / p{k)dk, 

M ^ /•■Bi 
T^Z^^ o-n(A)c?A, 

L ^ J-Bi 

M' ^ rB2 

T 



Vn / ujn{p)dp, 

n J-B2 

M" ^ fBs 



(7) 



orbital are given by, 

J pdk + '^n' J uJnd^ ~ X! y ^nd\ ^'^^ j Tndi 



L 2 

^£ = 1 
L 2 



j pdk j ujndjj.. (8) 

The energy is given by 

E f 

— = -2t J cos kp{k)dk, (9) 

and the magnetization by 

M, = gsS, + gtT„ (10) 

where Qs, Qt are Lande g factors that we know Qs = 2,gt = 1, which was ignored in previous 
hterature because it was not able to be related to spin-orbital model. 

III. THERMAL EQUILIBRIUM 

Yang and YangEl first introduced a definition of entropy in terms of distribution of roots 
and holes for Bethe-ansatz solvable systems. This method was recently carefully compared 
with the transfer matrix approach, the results of both approach are shown in agreement 
with each other0. In order to obtain thermal equilibrium at finite temperature, we should 
maximize the free energy Q = E + Ej — TS — AN where Ej = —^qHMz is the Zeemann 
energy due to the external magnetic field and A is the chemical potential. The entropy is 
given by 

l = Ji(p + Hp + p')-p\np-p'^ In p']dk 

+ E / ii^n + ^'J H^n + cr'n) ' In a„ - lu a^X 

+ E / [(^« + ^n) H^n + - In W„ - lu W^]ci/i 

+ ii^n + t!1) ln(r„ + rll) - r„ In r„ - In r„^]d^/. (11) 



We obtain the Helmholtz free energy from the condition = that 

T/ r 

F=-—Jln[l + e-^]dk, (12) 

where e should be solved from the following equations, 

e = [2t cosk + A+ ^HqH] - ^ / K„{\ - sin k) ln[l + e-'^-^^^jdX, 
On = 2fionH/T - J coskKnisink- A) ln[l + e~'^'''^]dk 

+ J2Ar,u f Kt{X-X') \n[l + e-'^^^'^]dX' -J2Bnit f Mfi - X) \n[l + e-^'^^^]dfi, 
It •' It 

Cn = -i2onH/T - Bnit [ Kt{X - /i) ln[l + e-'^^^^]dX 
It •' 

+ Anit f Ktifi - /i') ln[l + e-<^^^">]du' - Y Bnit f Kt{u - /i) ln[l + e-«'(^)]dz/, 
It •' It 

in = 2^IonH/T - J2 Bnit I Kt{fi - u) ln[l + e-^'^'^^jd/i 
It 

+ Ann I Kt{u - u') ln[l + e-^'^''''^]du' . (13) 
It 

In getting the above expressions, we have introduced notations: 

e^ = ^, e^" = ^, e^" = ^, e«" = ^. (14) 

P C^n Tn 

IV. THE CASE AT LOW TEMPERATURE 

At very low temperature, the density of roots undergo a slight modification from that of 
the ground state, and the right hand side of Eq. (1^) are approximately zero below Fermi 
surface. Then we have 



e = -— 2tcosk + A+ -jj^oH + / Ki{X ~ smk)e{X)dX, 

coskKi{smk- X)e{k)dk- / K2{X - X')9{X')dX' + / Ki{fi - X)C{fi)dfi, 

-Qo J-Bi J-B2 



Ki{X - fx)9iX)dX - K2ifx - p'XWl-i' + K,{iy - fx)i{iy)du, 

-Bi J-B2 J-B-i 

i = 2fioH/T + r K^{ii - z/)C(/i)rf/i - r K^iu - u')i{u')du'. (15) 

JBo. J-B'A 
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It is plausible to take Bi,B2,B3 = oo in our present case. By making a Fourier transform 
of the second through fourth equation of Eq. (p!5|) , we have 

e=2h + — cos fce-^l'"l-'"'^'"^e(fc)ciA; - e~2''l"'l^ + e'^'l^'IC, (16) 

271 J-Qo 

(= -h + e~"^''"^e - e-^^l^'IC + e-''l"'lf, 

where h = fioH6{w)/T. It is not difficult to derive an integral equation for e: 

e = -\2tcosk + A]/T+- t%osk'e{k')R^C-^^^^^^^^)dk (17) 

Tj J-kp ^ T] ^ 

where kp is the Fermi momentum and 

J_ sinh(nt.) 

In the absence of external field {H = 0), Eq.(|l^) is rigorous. When the external field 
is weak, it is also plausible because the changes on the integral interval of spin and orbit 
rapidity is of order of 1/L in the thermodynamic limit. Moreover, even if we consider the 
contribution of the densities outside the Fermi surface, it is of order e~^/^. Clearly the 
presence of external field has no explicit effect on the ratio of density of charge roots and 
holes at low temperature. We can also infer from Eg . ([T7|) and Eq. ([10|) there will be two 
critical fields H^, which can cause phase transition on the magnetization of the system. 

Once e is solved, the free energy at low temperature can be evaluated by Eq. (p!2|) . Though 
an explicit expression can not be obtained from Eq. (0) in general case, it is easier for 
numerical calculation. Define the density of roots at low temperature as p = po + Pi? where 
Pi is a slight deviation from the density of the ground state po0, we can obtain 

p = po 1 - e - - / ^ , . , TT^. (19) 

vr J + (sm k — A)^ 

Substituting e{k) of Eq. (p!7| ) into Eq. (p!2|), we obtain 



Clearly when T — the free energy is just the ground state energy except for the additional 
chemical potential term —AL and the e{k) just play a role of so called dressed energy. 

Strong coupling limit 

In strong coupling limit r] —>■ oo, we have 

e = -[2tcosk + A]/T, (21) 

Since we consider the problem at low concentration, the Fermi surface is much less than vr. 
Then we can rewrite Eq. (^T]) as 

e = [t/c2 _ /.]/T, (22) 
where k = 2t + A. And the distribution of the roots becomes 

P = ^l + e[*^'-«]/^' ^^^^ 



-kp 



with condition /^f^ p{k)dk = D, and when T = 0, it gives kp = \J Kq/I. We obtain 

(24) 



K = Kq 



2AkI 



where kq = tn^D'^ is the value of k at zero temperature. Obviously the system's energy 
becomes 



The specific heat and the entropy of the system are calculated as 



(25) 



Cv=S = ^, (26) 

which is Fermi-liquid like. 
Weak coupling limit 

In weak coupling r] —>■ 0, Eq.(|T^ can be solved exactly. In mathematics we have 

?;^o sinh(4?7) 4' 



Then the Eq. (|T8D is proportional to (5-function. The e{k) has the form 

e = -[8tcosfc + 4^]/T. (27) 

So we can repeat the steps in the section of strong couphng hmit. The specific heat and the 
entropy become 

Clearly it is still Fermi-liquid like. We can infer at low temperature the specific heat of 
the system is Fermi-liquid like regardless the value of coupling constant. The difference is 
the coefficient. 

Magnetic properties 

In the presence of external field, the densities of spin and orbit rapidities will change. 
We rewrite them as 

a = aoil-e'), 
LV = uJo{l - e'^), 
r = ro(l-e«). (29) 

where 6, (, ^ are all negative and proportional to H/T. Then the magnetization of the system 
has the form 

M,/L = ^Jpdk + J udfi -2 J ad\-2 J rdv. (30) 
Concerning M = 3N/4, M' = N/2, M" = N/4 at the ground state, we obtain 

MJL = 2 J aQC^dX - J cuoe^dfi + 2 J Toe^du. (31) 
and the susceptibility in the weak external field is 



X rp 



J aoe^dX - ^ / ^oe'-dfi + J roe^du. (32) 



However, if the external field is strong enough so that there are no spin and orbit flipped 
down. That is M = M' = M" = and the densities of A, /i, u are all zero. Then we have 
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p + p^ = p{l + e^) = l., (33) 



where e is given by 

1 r 



2t cos k + A + -fJ-oH 



2hqHJT = - y cos kK^{sm k - X)edk. (34) 

The He is the maximal value of {H^, H^} mentioned above. At large coupling constant the 
value of He scales like 1/r). Clearly, in the limit — > oo, a small value of H can flip all spin 
and orbit to upward. The density of charge roots becomes 

(35) 

which takes the value of l/27r at the ground state. The density of charge roots in strong 
coupling limit takes the same form as that in sufficient strong external field. 



V. SUMMARY AND ACKNOWLEDGMENTS 

In the above, we discussed onc-dimcnsional SU(4) Hubbard-hke Model in the thermo- 
dynamic limit, which is reliable at low temperature and low concentration. We studied the 
property of the thermal equilibrium and obtained the equilibrium equations. In the low 
temperature limit, the T-dependent density of roots and free energy, specific heat both in 
strong and weak coupling limit, as well as magnetic properties of the system are calculated. 
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